 &MODEL
 NLIST   =           1,
 NSTEP   =        2000,
 EPSI    =  1.0000000E-07,
 WGHT    =   1.000000    ,
 NOBS    =        1000,
 NVAXZ   =           4,
 NVAX    =           3,
 NORDER  =           1,
 NVAZ    =           1,
 NLAMZ   =           1,
 NLAMY   =           1,
 NLAMX   =           3,
 NSCALE  =           1,
 NEXPY   =           0,
 NRES    =           0,
 NCOR    =           0,
 NCOV    =           0,
 VARNAME = BC4HGAUPR,
 NVARIANT        =          24
 /
 &AUTOCOR
 RHO     = 2*0.000000000000000E+000  ,
 PI      =   1.00000000000000     ,  0.000000000000000E+000,
 NTPI    =           0,   538976288,
 MATFILE = Rdist100.in                                                                                                               
                                       
 /
 &HETERO
 NUMZ    =           4, 9*0,
 NUMLZ   =           1, 9*0,
 DELTA   = 10*0.000000000000000E+000  ,
 NTDEL   =           1, 9*0,
 VLZ     = 10*0.000000000000000E+000  ,
 NTLZ    =           1, 9*0
 /
 &YBOXCOX
 VLY     =  0.000000000000000E+000,
 NTLY    =           1,
 NUMLXY  =           0
 /
 &XBOXCOX
 NUMX    =           1,           2,           3, 47*0,
 NUMLX   =           1,           2,           3, 47*0,
 VLX     = 10*0.000000000000000E+000  ,
 NTLX    = 3*1, 7*0
 /
 &YXZDATA
 YXZFILE = TijL2.in                                                                        ,
 YXZFMT  = (5e15.7)                                                                        
 /

VARIABLE       MINIMUM         MEAN          MAXIMUM   NULL OBS.      SCALE
DEPENDENT
Tij         0.5899500D+04  0.2573965D+08  0.3419374D+09     0     0.1000000D-07
INDEPENDENT
Ti          0.1580900D+05  0.1860662D+06  0.3513510D+06     0     0.1000000D-04
Tj          0.6962000D+04  0.1851223D+06  0.4337820D+06     0     0.1000000D-04
Util_BC4    0.3868082D-07  0.7346229D-03  0.1399558D-01     0     0.1000000D+02
CONSTANT    0.1000000D+01  0.1000000D+01  0.1000000D+01     0     0.1000000D+01
HETEROSKEDASTICITY
TiTj        0.6827108D+09  0.3439253D+11  0.1093478D+12     0     0.1000000D-10


MRS of X(l) with respect to X(k) - at the sample means
===================================================================================================================================
        X(k) Ti          Tj          Util_BC4    CONSTANT  
X(l)
===================================================================================================================================

Ti           0.1000D+01  0.1229D+01  0.2776D+09  0.3413D+06

Tj           0.8137D+00  0.1000D+01  0.2259D+09  0.2777D+06

Util_BC4     0.3602D-08  0.4426D-08  0.1000D+01  0.1229D-02

CONSTANT     0.2930D-05  0.3601D-05  0.8135D+03  0.1000D+01


INITIAL ESTIMATES

             LAM    PARAMETER      STD-ERROR        STUDENT-T    GRADIENT       DERIV.Y         ELAST.Y
                                                                                DERIV.E(Y)      ELAST.E(Y)


Ti            X1  0.5325132D+00  0.6598434D-01      8.0702964  0.2799840D-11  0.7366572D+02      0.5325132
                                                                              0.1169046D+00      0.5325132

Tj            X2  0.6510872D+00  0.5906770D-01     11.0227274  0.1252758D-11  0.9052803D+02      0.6510872
                                                                              0.1475469D+00      0.6510872

Util_BC4      X3  0.5837007D+00  0.1616134D-01     36.1171024 -0.3518445D-10  0.2045165D+11      0.5837007
                                                                              0.1703872D+03      0.5837007

CONSTANT      X0  0.9767329D+00  0.1192295D+00      8.1920401  0.7217138D-11  0.2514076D+08      0.9767329
                                                                              0.3300557D+00      0.9767329

VAR(W)            0.1263045D+01  0.5648508D-01     22.3606798  0.0000000D+00

INITIAL VALUES
 
RHO1           0.0000000D+00
PI1            0.1000000D+01
 
TiTj       Z1  0.0000000D+00
LAM(Z1)        0.0000000D+00
 
LAM(Y)         0.0000000D+00
LAM(X1)        0.0000000D+00
LAM(X2)        0.0000000D+00
LAM(X3)        0.0000000D+00


LOG-LIKELIHOOD =  0.9344055D+03

NUMBER OF OBSERVATIONS              =           1000

PEARSON-R2(E) (unadjusted)          =         0.2435
PEARSON-R2(E) (adjusted for D.F.)   =         0.2412

PSEUDO-R2(L)  (unadjusted)          =         0.5780
PSEUDO-R2(L)  (adjusted for D.F.)   =         0.5767

PROBABILITY of Y to be at the limit =         0.0000
MEAN VALUE of E(Y)                  =  0.3379181D+00

MEAN of residual W                  =  0.9464949D-14
VARIANCE of residual W              =  0.1263045D+01
95% CONFIDENCE INTERVAL for residual W
Lower Bound = -0.2202751D+01   Upper Bound =  0.2202751D+01




initialization cards for control routine follow - -  id card - no. parameters , no. steps , step size 
minimization tolerance - starting values of parameters - relative uncertainties of parameters (error estimates) -  
and starting direction  (relative parameter increments for 1st step)
identification line
         7      2000   0.1E+01   0.1E-06
  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00
  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01  0.10000E+01
  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00  0.00000E+00

first entry to fcn , parameters of current step follow - - f,x(i)
  -0.9344054864D+03  0.00000D+00  0.00000D+00  0.00000D+00  0.00000D+00  0.00000D+00  0.00000D+00  0.00000D+00

new absolute minimum , parameters for fcn entry     5  follow - - f,x(i)
  -0.9453573924D+03  0.36907D-03 -0.20583D-01  0.00000D+00  0.63531D-01 -0.11152D-02  0.55656D-02 -0.20583D-01

new absolute minimum , parameters for fcn entry     8  follow - - f,x(i)
  -0.9511134348D+03 -0.32302D-03 -0.12195D+00 -0.28508D-02  0.83243D-01 -0.66533D-02  0.34873D-01  0.23977D-01

new absolute minimum , parameters for fcn entry    11  follow - - f,x(i)
  -0.9619196803D+03 -0.41515D-03 -0.36111D+00 -0.17324D-01  0.58543D-01 -0.20913D-01  0.11509D+00 -0.17241D-01

new absolute minimum , parameters for fcn entry    15  follow - - f,x(i)
  -0.9773634684D+03 -0.14084D-01 -0.44444D+00  0.37632D+00  0.42114D-01 -0.17525D+00  0.12554D+01 -0.41818D-01

new absolute minimum , parameters for fcn entry    18  follow - - f,x(i)
  -0.9800719720D+03 -0.18253D-01 -0.40781D+00  0.13345D+00  0.52703D-01 -0.23959D+00  0.15220D+01 -0.27054D-01

new absolute minimum , parameters for fcn entry    21  follow - - f,x(i)
  -0.9803403227D+03 -0.10149D+00 -0.41565D+00  0.14924D+00  0.54069D-01 -0.43663D+00  0.15191D+01 -0.29073D-01

new absolute minimum , parameters for fcn entry    25  follow - - f,x(i)
  -0.9814787574D+03 -0.77823D+00 -0.50011D+00  0.24324D+00  0.68871D-01 -0.35004D+00  0.16465D+01 -0.40575D-01

new absolute minimum , parameters for fcn entry    28  follow - - f,x(i)
  -0.9827911376D+03 -0.60334D+00 -0.68058D+00  0.50996D+00  0.75589D-01 -0.27287D+00  0.17886D+01 -0.26188D-01

new absolute minimum , parameters for fcn entry    32  follow - - f,x(i)
  -0.9852381772D+03 -0.22529D+01 -0.12868D+01  0.87070D+00  0.60601D-01 -0.44780D+00  0.17419D+01 -0.40470D-01

new absolute minimum , parameters for fcn entry    35  follow - - f,x(i)
  -0.9861824977D+03 -0.26369D+01 -0.15655D+01  0.10674D+01  0.73340D-01 -0.22829D+00  0.16615D+01 -0.29223D-01

new absolute minimum , parameters for fcn entry    38  follow - - f,x(i)
  -0.9866377223D+03 -0.37315D+01 -0.19800D+01  0.12173D+01  0.76664D-01 -0.33726D+00  0.15203D+01 -0.31033D-01

new absolute minimum , parameters for fcn entry    41  follow - - f,x(i)
  -0.9869569866D+03 -0.48323D+01 -0.24895D+01  0.15230D+01  0.75995D-01 -0.33267D+00  0.16153D+01 -0.33450D-01

new absolute minimum , parameters for fcn entry    44  follow - - f,x(i)
  -0.9871169478D+03 -0.54297D+01 -0.27887D+01  0.16051D+01  0.73582D-01 -0.31366D+00  0.16672D+01 -0.31378D-01

new absolute minimum , parameters for fcn entry    47  follow - - f,x(i)
  -0.9872573347D+03 -0.69426D+01 -0.33727D+01  0.18185D+01  0.76999D-01 -0.29825D+00  0.16176D+01 -0.35072D-01

new absolute minimum , parameters for fcn entry    49  follow - - f,x(i)
  -0.9872914704D+03 -0.76069D+01 -0.36699D+01  0.19146D+01  0.77683D-01 -0.33762D+00  0.16242D+01 -0.33896D-01

new absolute minimum , parameters for fcn entry    51  follow - - f,x(i)
  -0.9872958219D+03 -0.79512D+01 -0.38079D+01  0.19592D+01  0.77514D-01 -0.32709D+00  0.16259D+01 -0.34723D-01

new absolute minimum , parameters for fcn entry    53  follow - - f,x(i)
  -0.9872959398D+03 -0.79868D+01 -0.38222D+01  0.19624D+01  0.77664D-01 -0.32846D+00  0.16253D+01 -0.34681D-01

new absolute minimum , parameters for fcn entry    55  follow - - f,x(i)
  -0.9872959418D+03 -0.79938D+01 -0.38252D+01  0.19636D+01  0.77675D-01 -0.32857D+00  0.16253D+01 -0.34681D-01

new absolute minimum , parameters for fcn entry    56  follow - - f,x(i)
  -0.9872959418D+03 -0.79942D+01 -0.38254D+01  0.19637D+01  0.77674D-01 -0.32857D+00  0.16254D+01 -0.34682D-01

new absolute minimum , parameters for fcn entry    57  follow - - f,x(i)
  -0.9872959418D+03 -0.79943D+01 -0.38254D+01  0.19637D+01  0.77674D-01 -0.32857D+00  0.16254D+01 -0.34682D-01

new absolute minimum , parameters for fcn entry    58  follow - - f,x(i)
  -0.9872959418D+03 -0.79944D+01 -0.38255D+01  0.19637D+01  0.77674D-01 -0.32857D+00  0.16254D+01 -0.34683D-01


the function (fcn) has not improved by    0.10000D-06 after 3 steps which include a minimum - - -


therefore the run has been terminated
  -0.9872959418D+03 -0.79944D+01 -0.38255D+01  0.19637D+01  0.77674D-01 -0.32857D+00  0.16254D+01 -0.34683D-01


all numbers of last entry refer to final pt. of run


the lowest value of the function (fcn) was  -0.9872959418D+03 for entry   58 - - x(i) follows
 -0.79944E+01 -0.38255E+01  0.19637E+01  0.77674E-01 -0.32857E+00  0.16254E+01 -0.34683E-01



run statistics follow - -

 nfcn=   58 nst=   21


MRS of X(l) with respect to X(k) - at the sample means
===================================================================================================================================
        X(k) Ti          Tj          Util_BC4    CONSTANT  
X(l)
===================================================================================================================================

Ti           0.1000D+01  0.2036D+01  0.3179D+09 -0.2421D+08

Tj           0.4912D+00  0.1000D+01  0.1562D+09 -0.1189D+08

Util_BC4     0.3146D-08  0.6404D-08  0.1000D+01 -0.7617D-01

CONSTANT    -0.4130D-07 -0.8407D-07 -0.1313D+02  0.1000D+01


FINAL ESTIMATES

             LAM    PARAMETER      STD-ERROR        STUDENT-T      STUDENT-T    GRADIENT       DERIV.Y         ELAST.Y
                                   (CONDIT.)        (CONDIT.)      (CONDIT.)                   DERIV.E(Y)      ELAST.E(Y)
                                                    (PAR = 0)      (PAR = 1)

Ti            X1  0.8081107D+02  0.3293458D+03      0.2453684                -0.1125704D-08  0.5512397D+02      0.3984790
                                 0.1044874D+02      7.7340464                                0.8053409D+02      0.4245874

Tj            X2  0.8337020D-08  0.2854515D-07      0.2920643                -0.1396432D+00  0.1122142D+03      0.8070567
                                 0.6474783D-09     12.8761375                                0.1100343D+03      0.9053570

Util_BC4      X3  0.1465531D+01  0.5232246D+00      2.8009596                 0.4149658D-08  0.1752334D+11      0.5001251
                                 0.3225389D-01     45.4373424                                0.1741245D+12      0.5456669

CONSTANT      X0 -0.1951711D+03  0.7436661D+03     -0.2624445                -0.4361939D-09 -0.1334778D+10    -51.8568692
                                 0.3109953D+02     -6.2756923                               -0.1315232D+10    -52.3026395

VAR(W)            0.1878310D+02  0.1101520D+02      1.7051986                 0.0000000D+00
                                 0.1192805D+01     15.7470030
 
RHO1      (E)    -0.9999998D+00  0.6598055D+01     -0.1515598                -0.3807217D+00
                                 0.5962964D+01     -0.1677018
PI1       (F)     0.1000000D+01
 
TiTj      (E) Z1 -0.9596737D-21  0.1210415D-19     -0.0792847                -0.8681763D+16  0.2263949D-04      0.0302502
                                 0.1388170D-21     -6.9132295                               -0.1287221D-03     -0.2160342
LAM(Z1)   (E)     0.1963687D+01  0.5118953D+00      3.8361109      1.8825865  0.2056377D-03
 
LAM(Y)    (E)     0.7767373D-01  0.1847526D-01      4.2042020    -49.9222314 -0.1006926D-03
LAM(X1)   (E)    -0.3285703D+00  0.3528576D+00     -0.9311697     -3.7651741  0.3365654D-06
LAM(X2)   (E)     0.1625357D+01  0.2779381D+00      5.8479117      2.2499878  0.2884792D-05
LAM(X3)   (E)    -0.3468252D-01  0.1200557D-01     -2.8888696    -86.1835569 -0.1147465D-04


LOG-LIKELIHOOD = -0.1743338D+05

NUMBER OF OBSERVATIONS              =           1000

PEARSON-R2(E) (unadjusted)          =         0.3182
PEARSON-R2(E) (adjusted for D.F.)   =         0.3113

PSEUDO-R2(L)  (unadjusted)          =         0.9554
PSEUDO-R2(L)  (adjusted for D.F.)   =         0.9549

PROBABILITY of Y to be at the limit =         0.0000
MEAN VALUE of E(Y)                  =  0.2720163D+08

MEAN of residual W                  = -0.6715935D-02
VARIANCE of residual W              =  0.1878310D+02
95% CONFIDENCE INTERVAL for residual W
Lower Bound = -0.8501252D+01   Upper Bound =  0.8487820D+01


CPU TIME: 26 min 39 sec
